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We propose a Theoretical model accounting for the recently observed reduced magnetic moment 
of Gadolinium in fullerenes. While this reduction has been observed also for other trivalent rare- 
hearth atoms (Dy 3+ , Er 3+ , Ho 3+ ) in fullerenes and can be ascribed to crystal field effects, the 
explanation of this phenomena for Gd 3+ is not straightforward due to the sphericity of its ground 
state ( S= 7/2 , L= ). In our model the momentum lowering is the result of a subtle interplay 
between hybridisation and spin-orbit interaction. 



INTRODUCTION 

Endohcdral metallofullerenes M@Cs2 are novel materi- 
als that have attracted a wide interest in physics, chem- 
istry but also in material or biological sciences for the 
large variety of promising applications of their peculiar 
properties 0, 0, 0] In endohedral metallofullerenes, a 
positively charged core metal is off-center in a negatively 
charged strong carbon cage, resulting in strong metal- 
cage interaction and intrafullerene charge transfer from 
the metal to the cage. 0, H, 

The magnetism of these systems is mainly due to the 
spin of the entrapped metals. In a series of average mag- 
netization mesurements a paramagnetic behaviour has 
been observed la, III llfj , with negative Weiss temper- 
atures. The negative Weiss temperature indicates the 
presence of a weak antiferromagnetic interaction betwen 
the cage and the metal, and between neighbooring cages, 
but for heavy rare-earths (RE) endofullerenesflfj ferro- 
magnetic coupling has been mentioned in the sub-20K 
range. In the case of heavy RE these experiments gave 
a number of magnetons per encaged ion that is smaller 
than for the free ion. This result has been phenomenogi- 
cally ascribed to the cage crystal field interaction for high 
L ions and, for the L = Gd case, to the antiferommag- 
netic interaction between the ion and the cage. 

A recent work[ll] has used x-ray magnetic circular 
dichroism and (XMCD) x-ray absorption spectroscopy 
(XAS) to characterize local magnetic properties of heavy 
RE metallo centers, using the M4 and M5 resonances 
(3d — > 4/ transitions). The absorption spectra of this 
work were very well fitted assuming trivalent ions (4/™ 
electronic structure with n = 7 for Gd, n = 9 for Dy...), 
while XMCD confirmed that there is a strong reduction 
of the measured ion magnetisation compared to the free 
ion case. 

For L / ions the reduction was reproduced by a 
model Hamiltonian were a weak crystal field prevents the 
ion total angular moment J to align completely along the 
magnetic field. The case L — of trivalent Gadolinium 
was more difficult and hybridisation model gave not a 
satisfactory explanation. 



Infact, althought hybridisation gives antiferromagnetic 
coupling with the cage and accounts ( in the Gd case) for 
a 14% reduction of the average moment ( Gd + cages), 
it cannot explain the reduction of the moment localized 
on Gd ion (see next section). 

What we will show in the present paper is that a com- 
bined action of hybridisation and spin-orbit interaction 
can have a dramatic effect on the observed magnetic mo- 
ment. 

This effect is not trivial and is significant only in a 
restraint parameters region that was not discovered in 
the previous numerical study [llj. In the present paper 
we will give a complete analytical discussion of this ef- 
fect. In next section we will introduce a simple model 
with an anisotropic hybridisation where the Gd ground 
state is basically 4/ 7 with a small 4/ 8 component due 
to a backdonation from the cage, we will then consider 
spin-orbit interactions at the first order and show the 
dramatic changes in magnetisation. Analytical formula 
will be compared to exact numerical solutions. Finally 
we will discuss the possible experimental manifestations 
of the studied phenomena. 

HYBRIDISATION MODEL AND SPIN-ORBIT 

Back donation concerns a cage unpaired electron. As 
the Re ion is off center we consider that the hopping is 
non isotropic and choose to restrain the transfer of the 
cage backdonated electron to the 4/ orbital that is closest 
to the cage, the m z = one, where the z-axis is parallel 
to the off center displacement. 

The Gd ground state has S = 7/2, L — 0. Adding one 
more electron to this state, one can acces only the quan- 
tum numbers S = 3, L = 3 of the 4/ 8 configuration[12j. 

The energy difference between the 4/ 8 level and 4/ 7 
ground state will be named, without SO coupling, AE. 
This is a positive scalar quantity and must be large in 
magnitude, compared to the hopping strenght t, because 
the fractional backdonation has been observed to be very 
small. 

We consider an effective interaction term proportional 
to t 2 . This interaction transforms the state (S z , s) (where 
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S z is the spin z-component of the 4/ 7 shell and s the cage 
unpaired electron one) into the state [S z + 2s, — s), and 
the other way round. 

We can therefore restrict the effective Hamiltonian to 
the two-dimensional space spanned by (S z , s) and (S z + 
2s, — s ). The elements of the effective interactions are : 
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H(S Z , s —> S z + s — s ,s ) = 
{S z + s - s'\c 0s >\r))(r]\c£ s \S z 
E v — E 



(1) 



where Eq is the 4/ 7 ground state energy, c+ s and c ms are 
creation/annihilation operators for a 4/ electrons with 
m z = m and s z = s. The sum runs over the states 7/ of 
the 4/ 8 configuration. 

The case of Gd is quite simple because the only 
4/ 8 state having a non-zero parentage coefficient with 
|4/ 7 S* = 7/2L = 0) ground state is |4/ 8 S*' = 3, L = 3). 

The results of operating on the ground state on 4/ 7 can 
be expressed in terms of the parentage coefficients and 
angular rccoupling factors using well known formula |l3j|. 
for our specific case: 

c+ \4fSS z M z ) = -V8G S S ' L L ' x 
(SS z l/2s; S'(S Z + a)) x \4f 8 SL(S z + s)M z ) (2) 

Putting this formula into the Hamiltonian can be 
written as : 
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where vs z is the versor : 



v s , = {(7/2 + S,) l ' 2 ,-(7/2-S,+ 



(3) 



(4) 



For positive AE the ground state of H is v$ z , it cor- 
responds to antiferromagnetic alinement ( total angular 
moment J = 3), and has energy — ^gf ■ 

The state perpendicular to vs z has energy zero and 
corresponding total angular moment J = 4. There is no 
energy dependency on S z as it could have been expected 
on the basis of rotational invariance in the spin space, as 
long as spin-orbit interaction is not included. The local 
moment of Gd, in the antiferromagnetic ground state, 
can be almost fully aligned along the magnetic field. On 
the basis of equation Q one should observe, at saturation 



(S z ) = (7/2 x 7 + 5/2)/8 = 3.375 



(5) 



corresponding to a 3.6% reduction which is very far from 
the 20% observed one. The antiferromagnetic metal-cage 
coupling cannot explain the moment reduction observed 
with XAS, XMCD techniques. 

Therefore we add the spin-orbit (SO) coupling to the 
picture which breaks rotational invariance in the spin 



space. At first order SO splits the energies of the S — 3 
L = 3 state according to the total J , and affects the 
denominator involved in equation As a result the 
equation for H, that is given by equation for the zero- 
SO case, must be rewritten in this form : 



Hso(S z )=F(S z )v Sz ®v Sz 



where 



AE + E s f 



(6) 
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where S z = S z — 1/2. From a formal point of view 
SO gives at first order also another contribution beside 
affecting the propagator denominator. In fact the 4/ 7 
ground state, considering SO interaction, is not a pure 
S = 7/2 , L = state, but it has also a small component 
of the state S = 5/2, L = 1, whose amplitude is first 
order in SO strength. It is this perturbed ground state 
that should be considered in equation Q. However for 
symmetry reasons, in the framework of our model where 
backdonation affects only the m z = orbital, the con- 
tribution at first order in SO coming from this S = 5/2, 
L = 1 component is zero. 

At the moment we will therefore restrain the discussion 
to the above equation (JJJ) , that contains all the important 
physics of the studied phenomena. 

The S z dependency in equation (JJJ is given by the J 
dependecy of EjP. If the denominator in equation (7J is 
constant, one can factor the term : 
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and obtain again equation ©. For the 4/ 8 , S = 3, L = 



3 state the energy E^P is 



w so _ /■ 



24- J (J + f) 
12 



(9) 



where £ n z is the strength of SO interaction. This equa- 
tion can be obtained in a simple way : the state obtained 
putting seven spin-up electrons in the 4f shell, plus one 
spin-down electron in the m z = 3 state, has S — 3, L = 3, 
J = 6. 

It is very easy to calculate for this mono-determinantal 
state the expectation value of the S.L scalar operator, it 
is —1.5 ( m z — 3 of the eighth electron times its spin ). 
The expectation values for the other J can be calculated 
observing that the expectation value of a scalar product 
of two L = 1 tensors must be proportional to the 6J factor 
for the angular moments sextet (3,3, J ,1,1,0). By a quick 
glance at sixJ tables, equation JJJJ is readily obtained. 

The energy correction Ey, considering £„; = 
0.1975eV^4| has the negative value of about — 0.3eV for 
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J = 6 and a positive value of about 0.4eV for J =0. 
These values have to be compared to AE. Taking AE of 
the order of leV the effect of EjP is not negligeable and 
the dependency on S z is given mainly by the term with 
the smaller denominator, the J = 6 term. The Wigner 
symbol 

( 3 3 6 , ) (10) 

1 S z -s z ' y ' 

has the bigger value for the smaller S zl as one can under- 
stand classically considering that to get J = 6 one has 
to align S and L . 

Therefore the ground state has the smallest S z and, 
unless the polarising magnetic field is sufficiently strong, 
the observed local magnetic moment will be zero. 

The J = 6 term is the only one to consider for AE 
approaching the value of 0.3 eV, because its denominator 
in equation (JJJ) tends to zero. But a small denominator 
means strong hybridisation, while hybridisation is weak 
because the encaged Gd is an almost pure 4/ 7 configu- 
ration. 

One should therefore consider the region where AE + 
EjP is big compared to the hopping strenght. For AE 
going to infinity the studied effect cancels out( see equa- 
tion (£§). 

So we are going to study a region where the effect re- 
sults from a imperfect cancellation of the different terms 
involved in the sum of equation Q). 

As a cancellation is involved one has to be very precise 
evaluating each single term in the sum. Therefore we 
devote a particular attention to the exact values of EjP. 

These values could be obtained at second order using 
Racah formalism and summing contributions from all the 
4/ 8 states accessible operating with the spin-orbit inter- 
action on the 4/ 8 , S — 3, L = 3 state. However, for the 
scope of this paper, which is to clarify the consequences 
of equation J7J, it is sufficient to plug in the sum the 
energies obtained by exact diagonalisation of Gd 2+ ion. 

We show in table (1) the comparison of EjP energies, 
calculated at first order by equation JHJ, compared with 
the exact numerical result. 

The numerical value of EjP is obtained calulating nu- 
merically the energy eigenvalues of the 4/ 8 Hamiltonian, 
and subtracting the ground state energy of the 4/ 7 one, 
where SO interaction is accounted for in both Hamilto- 
nians. In the numerical calculation F2,F^,Fq are taken 
from Thole 0. The parameter Fq is already contained 
in AE, so we take Fq equal to zero in numerical calcula- 
tions. 

In figure (1) we show the energies as a function of AE 
for the antiferromagnetic eigenstates of H so (S z ) ( equa- 
tion © ) for S z between 1 (smallest energy ) and 3 (high- 
est energy), with a dashed line, compared to the exact 
numerical solution ( solid line) . The ground state energy 
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-0.295 


-0.098 


0.065 


0.197 


0.295 


0.361 


0.394 




-0.302 


-0.066 


0.076 


0.189 


0.268 


0.321 


0.346 



TABLE I: Dependancy of 4/ 8 ground states as a function of 
total moment J. First order formula ( first line) is compared 
to numerical results ( second line ) obtained using parameters 
from reference [T^| . Units are eV. 
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FIG. 1: Energies as a function of AE for S z between 1 (small- 
est energy ) and 3 (highest energy). The ground state for 
S z — is taken as origin of the energy scale. The dashed line 
is equation (JSJ with 4/ 8 energies energies at first order in SO 
(left) and calculated numerically for an isolated ion (right). 
Comparison is made with the exact numerical solution ( solid 
line) . 



for S z — has been subtracted (it is taken as origin of the 
energy scale). In the left panel the 4/ 8 energies entering 
equation© are calculated at first first order in SO, while 
in the right panel numerical 4/ 8 energies for an isolated 
ion are used. The parameters used in the calculation are 
t = 0.05el^ and AE between 0.4 and 1. One can observe 
that the simple formula [S] is an excellent approximation 
when exact energies are considered in the denominator. 

The energy splitting has to be compared with the mag- 
netic field strength. Considering a typical XAS-XMCD 
experimental case 0] with a 7 Tesla field, the energy 
gain got aligning about 7 Bohr magnetons from perpen- 
dicular to parallel direction with the field is 0.4meV. This 
energy is of the same order of magnitude of, or lower, 
than the splitting caused by hybridisation plus spin-orbit. 

This effect can therefore, depending on t and AE, pre- 
vail on the magnetic polarizing field and suppress par- 
tially, or completely, the magnetization. 
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FIG. 2: Magnetisation 1.8,3,6 and 8K temperature for the 
parameters choice AE — leV and t — 0.2eV. 



tal behaviour is reproduced. 

The above discussion leaves the problem still open. 
First of all further investigation is needed to better eval- 
uate the real values of AE and t, that in this work we 
have chosen arbitrarily with the only criteria of giving a 
numerical example based on conservative values ( small 
t and non negligeable AE). 

Second, a comparison with a complete set of data using 
a realistic modelisation of the sample should be done. 

However we can conclude that the effect is important 
even for very small hybridization, and therefore cannot 
be ignored if one wants really to understand the magnetic 
properties of encaged RE. 

I am grateful to the people of the ID08 beamline at 
ESRF, in particular Nick Brookes and Celine De Nadai 
for introducing me to this subject and motivating this 
analysis, and for the very fruitful discussions. 



DISCUSSION AND CONCLUSIONS 

We have shown in the previous section that a very 
small anisothropic hybridisation ( t = 0.05 eV ) can 
give, for weak magnetic fields, a complete suppression 
of the magnetisation along the encaged metal displace- 
ment axis. At zero temperature the magnetization would 
be a discontinuous function of the magnetizing field. For 
a polarizing field perpendicular to the displacement axis 
the magnetization curve would be instead continuous. 

In a real system one should take into account temper- 
ature and disorder. Temperature effect would tend to 
smear discontinuities. 

Disorder may be of different kinds. One kind of dis- 
order are fluctuations in AE due to inhomogeneities of 
cage environments. As AE affects the propagator de- 
nominator of equation (JZj) fluctuations might influence 
greatly the experimental result: discontinuities could be 
smeared out because the moment of cages having lower 
AE is depressed more than that of higher AE ones. 

Disorder of the displacement axis direction in the sam- 
ple would have a similar effect. 

These consideration could explain why experiments 
show magnetization curves that are continuous and sat- 
urate at reduced values. 

As an example we calculate magnetisation curves at 
different temperature in the case of random orientation 
of the displacement axis. We consider the parameters 
AE = leV and t — Q.2eV. The magnetisation is shown 
in figure 2 for 1.8,3,6 and 8K. These data have to be 
compared with figure 8 of reference|lO|. The experimen- 
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